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Abstract 

We generalize Wonham’s theorem on solvability of algebraic operator Riccati equations to 
Banach spaces, namely there is a unique stabilizing solution to A* P+PA—PBB*P+C*C = 0 
when (A, B) is exponentially stabilizable and (C,A) is exponentially detectable. The proof 
is based on a new approach that treats the linear part of the equation as the generator of a 
positive semigroup on the space of symmetric operators from a Banach space to its dual, and 
the quadratic part as an order concave map. A direct analog of global Newton’s iteration for 
concave functions is then used to approximate the solution, the approximations converge in 
the strong operator topology, and the convergence is monotone. The linearized equations are 
the well-known Lyapunov equations of the form A*P + PA = — Q, and semigroup stability 
criterion in terms of them is also generalized. 
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1 Introduction 

In this paper we generalize Wonham’s classical theory of positive solvability for algebraic matrix 
Riccati equations to Banach spaces. It provides a constructive solution to the linear quadratic 
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optimal control problem in Banach spaces that reduces to solving a sequence of linear operator 
equations [33, 12.3]. A generalization to Hilbert spaces was given by Zabczyk in 1970s [35], but 
for Banach spaces no general theory appears to exist despite the widespread use of the operator 
Riccati equations in the optimal control theory. For a Banach space X the operators appearing in 
the Riccati equation are not symmetric positive definite operators on X, they map from X to X*, 
but there are suitable notions of symmetricity and positive definiteness for them. The iterative 
process used in our solution goes back to Klcinman, who used it for matrix Riccati equations m, 
Wonharn in [33J 12.3] gave more general conditions under which it works. It consists of solving a 
sequence of Lyapunov equations of the form A*P + PA = —Q, where A is the generator of a Cq 
semigroup on a Banach space X, and Q : X — > X* is symmetric and positive definite. 

Two major issues complicate the solution theory on Banach spaces. One is the absence of 
positive de fini te isomorphisms, like the identity operator on a Hilbert space, which allow non¬ 
degenerate approximations. These can be sidestepped however. The second issue is deeper. The 
Riccati equations involve operators on spaces of operators that lack regularity properties in the 
the strong and the weak operator topologies. It turns out that the suitable topology is the weak*, 
and explicitly or implicitly one has to work with Banach preduals to spaces of operators. For 
Hilbert spaces the predual can be identified with the space of trace class operators, but in the 
Banach case it is described as a projective tensor product. Tensor products of Banach spaces are 
rarely used and relatively little known in control theory applications, which helps explain why the 
operator Riccati equations are either considered in Hilbert spaces only [1, [36], or in some special 
cases [I3l [251126]J . 

We give a novel interpretation of infinite dimensional Riccati equations in terms of Lyapunov 
semigroups and their generators, see [8j for a similar approach with matrices. As pointed out in 
HB] and [19] the theory of Lyapunov equations implicitly relies on the fact that the left hand side is 
the generator of a semigroup T(t)P := T*(t)PT(t ) on the space of symmetric bounded operators 
B S (X,X*). We called it the Lyapunov semigroup, and one of its attractive properties is that it is 
positive, i.e. it preserves the positive definiteness of P. Unfortunately, this semigroup is not Co 
if A is unbounded, even when X is a Hilbert space, which obstructs its use for analytic reasons. 
However, in [T9] we proved that it is always adjoint to a C 0 semigroup on a predual to B S (X,X*), 
i.e. T{t) is always a Cq semigroup P, 3.1]. This property allowed us to overcome the analytic 
obstructions, but the price for using it is the appearence of tensor products as preduals. We should 
note that the original motivation for studying Cq semigroups came from automorphism flows on 
von Neumann algebras, a particular case of Lyapunov semigroups when T(t ) is unitary. 

In this paper we extend an explicit description of the Lyapunov semigroup generators D2] to 
non-reflexive Banach spaces, see Theorem [51 and use it to solve operator Riccati equations. In 
addition, we also generalize a stability criterion in terms of Lyapunov equations [19] to non-reflexive 
spaces, see Theorem [U the matrix case is due to Wonharn [33], 12.4], 

That Riccati equation can be written as AP + < h(T > ) = 0, where A generates a positive Cq 
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semigroup, is only half of the solution however. The other half comes from the fact that the 
quadratic part 4>(.P) is Gateaux differentiable and concave in the sense of partial order on the 
space of operators, a generalization of matrix concavity [3]. It explains why Kleinman’s method 
led to monotone convergence of approximations, as shown in [8] it amounts to an analog of Newton’s 
iteration on the space of matrices. 

Applications to control theory single out a particular class of solutions to the equation that 
translate into asymptotically stabilizing controls. It turns out that such stabilizing solutions P 
can also be characterized abstractly, for them the formal derivative A + $'(P) generates a stable 
semigroup. We found it more illuminating to be general and prove existence of stabilizing solutions 
for a class of general quasi-linear concave equations, Theorem [61 Aside from admitting a more 
streamlined argument, the result is of interest in its own right because it may apply to other 
situations, e.g. to reaction-diffusion equations with concave non-linear parts [22l 4.2], After 
everything is said and done our final result can be stated (almost) identically to Wonham’s [33], 
12.3]: if A — BK and A — LC generate exponentially stable semigroups for some I\ and L, the 
algebraic Riccati equation A*P+PA — PBB*P+C*C = 0 has a unique positive definite stabilizing 
solution P , see Theorem [2l 

The paper is organized as follows. Section [21 serves as preliminaries, we introduce the main 
concepts of the paper and state our main results. In particular, the linear quadratic optimal 
control problem in a Banach space X is recalled, and how it leads to the algebraic operator Riccati 
equation. Analytic difficulties related to its interpretation and solution are also discussed. Sections 
[3] and [4] lay the analytic groundwork for the proof which is completed in Section [5} In particular, 
Section [3] introduces projective tensor products, which are preduals to the spaces of bounded 
operators between Banach spaces, considers topologies in which the Lyapunov semigroups are 
continuous, characterizes their generators, investigates properties of positive definite symmetric 
operators, and continuity of quadratic maps on monotone sequences of operators. In Section [4] we 
prove an abstract result on existence of stabilizing solutions to quasilinear concave equations in 
ordered spaces, and in Section [5] we deduce from it our generalization of the Wonham’s theorem. 

2 Optimal control and operator Riccati equations 

This section provides motivation for the theory developed in the paper and states our main re¬ 
sults. We start by describing the optimal control problem that leads to algebraic operator Riccati 
equations, then introduce the main concepts needed to interpret and solve them, and conclude 
with discussion and examples. 

The linear quadratic regulator problem with infinite horizon can be described as follows. Let X 
be a real Banach space, called the state space, and A be the generator of a Co semigroup on X with 
domain T>a- Let U and V be real Hilbert spaces, called the control (input) and the observation 
(output) space respectively, and let B : U —* X and C : X —> V be bounded linear operators. 
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Consider a linear system written formally as 


x = Ax + Bu , x(0) = Xq 
y = Cx. 


Here u : [0, oo) —> U is a control input meant to steer the state x(t) so as to fulfill some objective. 
The state itself is not available for direct observation, only some reduction of it y{t) is. One typical 
control objective is optimal asymptotic stabilization, that is driving the state towards the origin 
asymptotically while minimizing a ’’cost”. A popular way to do so is to minimize a quadratic cost 
functional 

poo 

j[u,x o] := / (y,y)v + {u,u)udt ->• min, (2) 

Jo 

where (,)u, (, )v are the inner products on U and V respectively. By analogy to the Hilbert 
case [ 33l 12.1], [36] IV.4.1], one can show formally that the optimal control can be found in the 
feedback form u(t) = —B*Px(t), where P is a bounded symmetric positive definite operator that 
solves the algebraic operator Riccati equation: 


A* P + PA — PBB*P + C*C = 0. 


(3) 


Note that when the feedback control is applied, again formally, the state evolution equation be¬ 
comes x = (A — BB*P)x, so the control objective can be accomplished only if A— BB*P generates 
a stable semigroup. When X is a Banach space P has to be interpreted as an operator from X to 
X*, but the usual Hilbert notions generalize naturally to such operators. 

Definition 1 . Let B(X,Y) denote the Banach space of bounded linear operators from a Banach 
space X to a Banach space Y with the induced norm. Denoting (, ) the duality pairing between 
X and X* we call P G B(X,X*) symmetric if (Px,y) = (Py,x) for all x,y G X, or equivalently 
P*\ = P (where X is identified with a subspace of X**). We call P positive definite if (Px , x) > 0 

for all x G X. The subspace of symmetric operators is denoted B S (X. X*), and the cone of positive 
definite ones in it Bf(X,X*). Notations B S (X*,X) and Bf(X*,X) are defined similarly. 

Most authors dealing with operator Riccati equations assume X to be Hilbert (see however 
[231 [251 [26]). in which case X* can be identified with X, and P is a self-adjoint operator on A". 
However, in applications the state space is dictated by the problem, and is often not a Hilbert 
space, ft may seem natural to also allow the control and the observation spaces U and V to 
be non-Hilbert, but that leads to no real generalization as long as the cost functional remains 
quadratic. If one simply replaces the inner products in (J2D by positive definite quadratic forms we 
can use them to define inner products, which brings us back to the Hilbert spaces U and V. 
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As usual, we identify U and V with their duals, and set N := BB* £ Bf(X*,X) and Q := 
C*C £ Bf(X, X*) for brevity. Since A and A* are only densely defined some care is needed to 
define solutions to (|3]) . The traditional interpretation is to require that for all x, y £ T>a one has 

(Px, Ay) + (PAx, y) - ( Px , NPy) + {Qx, y) = 0. (4) 

In the Hilbert case this is called inner product Riccati equation [36] IV.4.2], so it is natural in 
general to call it the duality pairing Riccati equation. It is considered e.g. in [23] (25, 2SJ. We will 
use a different interpretation of (J3J) which treats P t—)■ A*P + PA as the generator of a positive 
semigroup. 

Definition 2. Let T(t) be a Co semigroup on X with the generator A. The semigroup T{t)P : = 
T*(t) PT(t) on J3(X,X*) and its restriction to B S (X,X*) will be called the Lyapunov semigroup 
ofT(t), and its generator will be denoted La and called the Lyapunov generator. 

The semigroup property of T(t) is straightforward to verify as is its positivity, indeed ( (T(t)P) x, x) = 
(PT(t)x,T(t)x) > 0 for positive definite P. But it is not obvious in general that it is continuous 
in some useful sense, for instance it is never Co for unbounded A. However, it was shown in ra 
that Lyapunov semigroups are adjoint to Co semigroups on preduals to B(X,X*) and B S (X,X*), 
in other words they are always Cq semigroups |5] 3.1]. 

This allows us to make sense of the generator La and use its standard properties. Moreover, it 
turns out, see Theorem [3] that the domain of La consists of operators P that satisfy P{T>a) Q Pa* 
and have ||A*P + PA\\ < oo. Then LaP is the extension of the bounded operator A*P + PA from 
T>a to the entire space. Therefore, we can interpret equation (j3j) as AP + <E>(P) = 0 with A = La 
and $(P) = — PNP + Q, and view solutions as elements of T>a that satisfy it literally. It also 
follows from Theorem [3] that solutions so defined are the same as solutions to the duality pairing 
Riccati equation. Our conditions for their existence are in terms of stabilizability and detectability, 
which for matrices go back to Wonham [331 12.6]. 

Definition 3. Let A be a generator of a Co semigroup on X and B e B(U,X), C G B(X,V). 
The pair (A, B ) is called exponentially stabilizahle if there exists I\ G B(X, U ) such that A — 
BI\ generates an exponentially stable C 0 semigroup, and the pair (C, A) is called exponentially 
detectable if there exists L £ B(V,X) such that A — LC generates an exponentially stable Cq 
semigroup. 

Note that disregarding continuity stabilizability is formally dual to detectability, i.e. ( A , B) is ex¬ 
ponentially stabilizable whenever (B*,A*) is exponentially detectable. For more on their meaning 
in control theory see Section [5} 

We will now state our theorems for Lyapunov and Riccati equations. In the following we abuse 
notation by identifying operators A*P + PA with their closures. 
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Theorem 1. Let X be a Banach space and T(t ) be a Co semigroup on it with the generator A. If 
the pair ( C , A) is exponentially detectable then the following conditions are equivalent: 

(i) A*P + PA = —C*C has a positive definite solution P £ Bf(X,X*) such that P(T>a) F Pa*', 

(ii) T(t ) is exponentially stable; 

(iii) The Lyapunov generator £ a has a bounded inverse on B S (X,X*) and — (Ca)~ 1 > 0. 

Recall that we defined ” positive definite” by a non-strict inequality, so for (7 = 0 the unique 
solution is P = 0, and the theorem is vacuously true (T(t) is exponentially stable by definition of 
(0, A) being exponentially detectable). The exponential detectability above can be replaced with a 
weaker condition of detectability in L 2 , see Section 0 This theorem was proved in [T9] for reflexive 
spaces, but the reflexivity assumption can be dropped using Theorem [3] 

Theorem 2. Let A be the generator of a Co semigroup on X, U,V be Hilbert spaces, and B £ 
B(U, X), C £ B(X, V). Suppose that ( A , B) is exponentially stabilizable and ((7, A) is exponentially 
detectable. Then 

(i) The algebraic operator Riccati equation 

A* P + PA — PBB*P + C*C = 0 (5) 

has a unique positive solution P such that P(Pa) Q Pa*, and ||7l*P + P7l|| < oo. This solution is 
stabilizing, i.e. A — BB*P generates an exponentially stable semigroup . 

(ii) Pick P 0 so that A — BB*P 0 generates an exponentially stable semigroup. This can always 
be done, e.g. by solving the Lyapunov equation 

(A - BI<)*P 0 + P 0 (A - BI<) = —C*C - K*K , 

with A — BK generating an exponentially stable semigroup. Then the solution P can be obtained 
as a monotone strong operator limit of solutions to linear Lyapunov equations 

(A - BB*P n )*P n+ 1 + P n+1 (A - BB*P n ) = —C*C - P n BB*P n . (6) 

Moreover, there is k > 0 such that ||P — P n +i|| < x|| P — P n \\ 2 (quadratic convergence). 

The approximations in part (ii) are obtained by formally applying Newton’s method to the equation 
F(P) := CaP + d>(P) = 0, and monotonicity is due to the fact that F(P) is order concave, that 
is F(aP + (1 — a)R) > aF(P) + (1 — a)F(R) for a £ [0,1]. The proof of convergence relies on 
the properties of Cq generators. In fact, we first prove an existence result in a general setting 
of order concave equations with linear parts generating positive Cq semigroups, see Theorem [6] 
Quadratic convergence, unfortunately, is not nearly as useful in infinite dimensional Banach spaces 
as in finite dimensional ones. The theorem only assures strong operator convergence of P n to P , 
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so || P — P n \\ may never get smaller than -. But only if that happens does quadratic convergence 
imply exponential rate of convergence by norm. 

Note that that B and C enter ([5]) only through N := BB* G Bf(X*,X) and Q := C*C G 
£>+ (X. X*). One could restate the theorem using only these operators because they can always be 
canonically factored. For Q for example define an inner product on Ran(Q) c X* by ( Qx,Qy) : = 
( Qx,y ). This is well defined, if Qx = Qx and Qy = Qy the value is the same since Q*\ x = Q- 
Denote by V the completion of Ran(Q) in the Hilbert norm, this V is called the reproducing kernel 
Hilbert space [13], [23] - Set Cx := Qx , to find the adjoint take z = Qy G Ran(Q) and compute 

(Cx, z) = ( Qx, Qy) = ( Qx, y) = (x, Qy) = (x, z). 

So C* acts on Ran(Q) as the natural inclusion to X*, and hence coincides with it on V since the 
range is dense. Thus, Q = C*C, and C is canonically recovered from Q. When X is a Hilbert space 
C can be identified with the positive square root Q^ , and V with Ran(Q 2 ). A similar construction 
factorizes N. To illustrate the theorem we give a couple of examples. 

Example 1: Any operator N G B+(X*, X) defines multiplication on B(X, X*) by (P, R) (->• PNR. 
Given also a Q G BJ~(X, X*) we can talk about taking a ’’positive square root”, i.e. solving 
PNP = Q for P G B+(X, X*). It is clear however that even for matrices this equation is not 
solvable for all pairs N and Q. 

Consider instead a regularized equation PNP + 2 aP = Q for some a > 0. Written as —2 aP — 
PNP + Q = 0 it is of the form ([5]) with A = — al , where / is the identity operator. Since this 
A generates an exponentially stable semigroup for any a > 0 the pair (A, B) is exponentially 
stabilizable for any B, and ((7, A) is exponentially detectable for any C. Theorem [5] now implies 
existence of a unique stabilizing solution. One can take Po := ^ l -Q as the initial guess for the 
Newton’s iteration. 

Note that for a = 0 Wonham’s theory does not apply in general even if PNP = Q is solvable. 
Existence of a stabilizing solution implies that — NP generates an exponentially stable semigroup, 
and therefore is invertible. Thus, a linear isomorphism of a Banach space would have to factor 
through a Hilbert space, a very special property. 

Example 2: In many applications the control and the observation spaces are finite dimensional 
[Tj . For asymptotic stabilization to be possible the ’’unstable part” of the generator A has to be 
finite-dimensional as well. To make this more precise assume that A is exponentially dichotomous, 
i.e. the space A" is a direct sum A_+A + of T(t) invariant subspaces such that the restrictions 
A_ := A\ x and A + := — A\ x generate exponentially stable semigroups [56]. It follows that 
the spectrum of A does not meet the imaginary axis, and if A is bounded this suffices for the 
exponential dichotomy. 

Assume that dim(A + ) < oo and X + C RanP. Denoting I + the identity operator on X + we 
can use a basis in X + to define K G B(X, U ) such that BK = 0 + (A + + I + ). Then A — BK = 
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A _ + (—1+) generates an exponentially stable semigroup, and {A, B) is exponentially stabilizable. 
Similarly, (C,A) is exponentially detectable if KerC C X_ since then we can define L £ B(V,X) 
such that LC = 0 + {A + + I + ). Thus, if the range and the kernel conditions are satisfied there is 
a unique stabilizing feedback control, a result widely used in Hilbert spaces. 

Note that we actually need something weaker than the exponential dichotomy in this case, as 
long as dim(X + ) is finite A + can be allowed to have purely imaginary or even positive real part 
spectrum. If dim( X + ) = oo however the range and the kernel conditions are not enough, and their 
replacements depend on conditions for factorization through Hilbert spaces [22J Ch.2], 


3 Topologies and continuity 


This section develops some analytic tools for studying the operator Riccati equations, and proving 
convergence of approximations to their solutions. First, we introduce projective tensor products 
that are preduals to spaces of operators. Then we consider implemented semigroups, that restrict 
to the Lyapunov semigroups, and tensor product semigroups, that they are adjoint to. Their 
continuity in some natural operator topologies is discussed, and the generators are described. To 
handle the nonlinear part of the equations we then look into properties of positive cones in spaces 
of symmetric operators, and continuity of quadratic maps on monotone sequences in them. 

Definition 4. Let X and Y be Banach spaces and X <g> Y be their algebraic tensor product. The 
duality pairing between B(X,Y*) and X ® Y is defined by ({x <8) y, P)) := ( Px,y) on monomials, 
and extended by linearity. Given p £ X ®Y its projective tensor norm is defined by 


:= inf'{yb \\x t 


P = Ei x % ® Vi, x i e X, yi £ Y} . 


The projective tensor product X <8^ Y is the completion of X <8> Y in this norm [28, 2.1 ]. 

The projective tensor product is convenient for several reasons. First, its elements can be described 
explicitly, they are of the form a i x i®Vi with bounded sequences x t , y t £ X , and a summable 
numerical sequence a* £ R (this is the Grothendieck representation theorem [3TJ III.6.4]). Second, 
since ||x £8 y\\ = ||x||||?/|| it is almost immediate that the pairing ((•, •)) extends from X <g) Y to 
X<S> n Y, giving us a Banach space isomorphism ~ B(X, H*). Moreover, the dual norm on 

B(X,Y*) is exactly the induced operator norm, so this isomorphism is isometric [28i 2.1]. Finally, 
the semigroups on X (g)?,. Y that the Lyapunov semigroups are adjoint to also admit an explicit 
description, see Definition 0 

Now let us recall some natural topologies on operator spaces. Given two Banach spaces X and 
Y consider the space B(X, Y*) of bounded operators from X to the Banach dual of Y. Aside from 
the weak* topology induced by duality with the predual X <8^ Y the following three topologies 
will be of use to us: 



• Weak operator* (wo*) with seminorms R \(Rx, y) \ for x E X, y E Y\ 

• Strong operator (so) with seminorms R i—>■ ||i?x|| for x E X\ 

• Ultraweak* (uw*) with seminorms R n- | {R x i-, Vi) I f° r x i G X,yi £ Y, and 
ESi W 2 < 00, ESi IM 2 < 00. 

The weak operator* topology is obviously weaker than the other two, and the same argument 
as for Hilbert spaces [9j 1.3.2] shows that it coincides with the ultraweak* topology on bounded 
subsets of B(X,Y*). Recall that given a locally convex space X and a separating subspace y of 
its algebraic dual, cr(X, y) denotes the weakest topology on X in which all functionals from y are 
continuous, ft is specified by the seminorms u i —> |(u, u)|, v E y. With this notation we have the 
following. 

Lemma 1. The weak operator* topology coincides with <j(B(X,Y*),X <E> Y), and the ultraweak* 
topology coincides with the weak* topology, a(B(X,Y*),X <8v Y ). 

Proof. The first claim is almost obvious from definitions since (Rx, y) = (( x®y , R)). For the second 
claim let p(R) | Ui) \ be an ultraweak* continuous seminorm, and set p := x i®Vi- 

Taking the projective tensor norm: 

OO OO OO OO 1 OO 1 

iipii : = = £mni - (5Zimi 2 ) 2 (5^ii^ii 2 ) 2 <(X) ' 

i —1 i =1 i =1 i =1 i =1 

Thus, p E X r Y and p(R) = \((p, R))\ is a a(B(X, Y*),X Y) seminorm. 

Conversely, let p E X Y. By the Grothendieck representation theorem p = a* Ui ® Vi 

with 11^11,11^11 < M < oo and |a,| < oo [HD III.6.4]. Set Xi := sign (a,) | a* | 2 Uj and yi := 
\ai\iVi, then ||ay|| 2 < oo, ET= i INI 2 < oo, and |((p, R))\ = \ Y^=i{R x ii y%) | is an ultraweak* 
seminorm. Thus the ultraweak* and the weak* topologies share the same continuous seminorms, 
and hence coincide. □ 

There are two closely related classes of semigroups on B(X, Y*) and X <g)„. Y induced by semi¬ 
groups of bounded operators on X and Y, see [ID 1.3.16], [SDJ 3.4] and [12]. 

Definition 5. Let X and Y be Banach spaces and T[t ) and S(t ) be Cq semigroups on X and 

Y respectively. The semigroup T(t)P := S*(t) PT(t) on B(X,Y*) will be called the semigroup 
implemented by T(t) and S(t). For a topology r on B(X, Y*) we define the r-generator of Tit ) by 
A T P := r-lim^o f (T(t)P — P ) on the domain Vjt where the limit exists. The semigroup defined 
by 7 *(t)(x <g) y) := S{t)x <g) T{t)y mid extended by linearity and continuity to X (gv Y is called the 
tensor product semigroup induced by S(t ) and T(t). 
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As the notation indicates T(t) is adjoint to 7*(t) [TH1, Thm.4], and therefore is a weak* continuous 
or Cq semigroup on B(X,Y*). By Lemma Q] it is also ultraweak* and hence weak operator* 
continuous. Since the last two topologies coincide on bounded subsets we have M uw = M wo , and 
this generator is weak* densely defined and weak* closed. In general, S*(t) may only be a Cq 
semigroup, but if Y is reflexive it is a Cq semigroup on Y* (5, 3.1.8], [33] IX.1], In that case T(t) 
is also strong operator continuous. Indeed, 

\\(T(t + h)P-T(t)P)x\\ 

< ||5*(t + h) || ||P|| || T(t + h)x - T(t)x || + || (S*(t + h)- S*(t))PT(t)x || -> 0 . 

h—>0 

Moreover, T{t) is locally uniformly bounded in the sense of [18], and locally bi-continuous in the 
sense of Kiihnemund [20[ Prop.3.16]. Either property implies that M so is strong operator densely 
defined, and strong operator sequentially continuous. 

As observed by Freeman ra, when Y is reflexive the strong operator generator coincides with 
the weak operator one, and can be described explicitly. We now prove a more general result that 
applies to non-reflexive spaces, and also describes the generator in terms of the duality pairing. 

Theorem 3. Let T(t ) be a semigroup on B(X,Y*) implemented by Cq semigroups T(t ) and S(t) 
on X and Y respectively, with generators A and B respectively. Define AP := B*P + PA (the bar 
stands for closure) on the domain 

V A ■= {P G B(X, Y*) | P(V A ) C XV, || B*P + PA\\ < oo} . 


Then, 

(i) P 6 V A if and only if there exists a bounded operator Q e B(X, Y*) such that 
(Px, By) + ( PAx , y) = (Qx, y) for all x G V A and y e V B , in which case AP = Q; 

(ii) A wo * = A; 

(iii) If additionally Y is reflexive then A so = A; 

Proof, (i) Assume that such Q exists. Let x G D A and z Px, then for any y G V B the definition 
of Q yields 

(A By) = -(PAx, y) + (Qx, y) = ((-PA + Q)x, y). 

By definition of the adjoint it follows that z = Px G V B * and B*(Px) = —PAx + Qx. Therefore, 
P(T>a ) Q Pb* and || B*P + XMj| = ||Q|| < oo. Thus, P G V A and AP = Q. 

Conversely, if P G V A then AP is bounded and AP = B*P + PA on V A . Taking x G V A and 
y G V B we have (( AP)x, y) = (Px, Gy) + (PAx, y), so Q := AP exists. 
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(ii) First we prove that A C n. wo *. Let P E V A and x E T>a , y G Y, then 

\{ (r(t + h)p - r(t)p)x, y) = \( (s*(t + h)pr(t + h)- s*(t)PT(tj) x , y ) 
= (s*{t + h)pt (V(f + h) - T(*))x, y) + Q (s*(t + h) - S*{t)) PT(t)x, y) 
= (^(t( t + h)~ T(tj)x, P*S(t + h)y } + (l (s*(t + h) - S*{t))PT(t)x, y } . 


Since x G also T(t)x G Pa, and in the first term \{T[t + h ) — T(t))x —> AT(t)x by norm. 
Moreover, P*S(t + h)y » P*S(t)y by norm since /S'(t) is a Co semigroup. Therefore, as h —> 0 

h —>-0 

the first term converges to ( T(t)Ax , P*S{t)y) = ( S*{t)PT(t)Ax , y). Since T(t)x G Pa and P G Pa 
we have PT{t)x G P_b*. Hence the second term converges to ( S*(t)B*PT(t)Ax,y). Summarizing 
we see that 

j t (T(t)Px,y) = (S*(t)PAT(t)x, y) + ( S*(t)B*PT(t)x,y) = (S*(t)(B*P + PA)T(t)x,y ). 
Integrating both sides from 0 to h we obtain, 

{(T(h)P-P)x,y)= f ( S*(t)(B*P + PA)T(t)x,y)dt. (7) 

Jo 

Since P G Pa the operator B*P + PA is bounded on Pa, and it has a unique bounded extension 
to X, which is its closure B*P + PA. The equality in (J7J) extends to all x G X if B*P + PH 
is replaced by P*P + PA. Dividing both sides by h and passing to limit as h —» 0 we arrive at 
{(A wo *P)x,y) = {(B*P + PA)x,y) for all x E X, y E Y. So Pa C V a ^, 0 * and A wo *P = AP on 
Pa, i.e. A C M wo *. 

Now we prove that M wo * C M. Take P eV a wo* cind consider th.G elemenljnry identity 

S-(h) -l)P= l(T(h)P ~P)~ S-(h)P^ (T(h) - /) . 

Applying both sides to x E V A and taking the weak limit in X as h —» 0 we obtain on the right 
(A wo *P):x — PAx. Therefore, Px E Vg* and B*Px = ( A wo *P)x — PAx. Hence P(Pa) Q Pb* 
and A wo *P = B*P + PA on V A ■ Since A wo *P is bounded so is B*P + PH, and since V A is dense 
A wo *P = B*P + PA everywhere. Thus, P G V A and H wo * = A on V A ^ 0 *. Combining the two 
inclusions we have H. wo * = A. 

(iii) When Y is reflexive the weak and the weak* topologies on it coincide, S*(t) is a Co 
semigroup, and its weak generator B* is the same as its strong generator by a theorem of Yosida 
p£h IX. 1], Hence we immediately have from (ii) that H so C AA°* C A. 


11 








For the converse inclusion note that x £ T>a still implies PT(t)x £ T>b* and (J7J) holds. But now 
S*(t)(B*P + PA)T(t)x is norm continuous, and (T(h)P — P)x = S*(t)(B*P + PAf)T(t)x dt is 

a Bochner integral. Reasoning as in (ii) we now get A C .4 s0 and A = M so . □ 

When Y is not reflexive T(t) may not be strong operator continuous, and T>a so may be strictly 
smaller than T>a- In other words, for some P £ T> 4 the difference quotient \{T{h)P — P) may not 
strong operator converge to B*P + PA. 

For the Lyapunov semigroups A = Y, T(t ) = S(t), and we restrict T{t) to the subspace of 
symmetric operators B S (X,X*), a predual to which is the subspace of symmetric tensors. 

Definition 6. The algebraic symmetric tensor product X®X is the linear span in X<S>X of tensors 
of the form x ® y + y ® x, where x, y £ X, and the projective symmetric tensor product X<g> n X is 
the closure of X®X in X X. 

The projective duality (X^Y)* ~ B(X, Y*) [28] 2.1], and the standard formulas for isomorphisms 
of subspaces and quotients [2914.8] show that also (Acgi^A)* ~ B s (X. X*). By inspection, Theorem 
[3] remains true for the restricted semigroups without change. It then follows from Theorem |3](i) 
that the duality pairing Riccati equation ([4]) is equivalent to the Lyapunov generator equation 
CaP + < L(P) = 0 in the sense that any solution to one is also a solution to the other. The following 
property of Lyapunov generators is immediately obvious from their explicit description in the 
theorem, and will be used in Section 0 

Corollary 1. Let X be a Banach space, A be a generator of a Co semigroup on X, and Ca be the 
corresponding Lyapunov generator. Then for any and G £ B(X) we have Ca+g = £a + C-g, in 
particular Vc A+G = Pc A ■ 

Our second concern is the continuity of the quadratic part of the Riccati equation. To this end 
we investigate relations between topology and order in B s (X, X*). When X = H is a Hilbert space 
B S (X, X*) = B S (H) is the familiar space of bounded self-adjoint operators with the cone of positive 
definite operators defining the partial order. As in the Hilbert case B S (X,X*) is isomorphic to 
the space of bounded quadratic (or bilinear symmetric) forms on X [28l 2.2], For the quadratic 
forms the partial order is just the pointwise order. However, in some respects the partial order on 
B S (X,X*) may behave quite differently than on B S (H). 

For one, unless A" is isomorphic to a Hilbert space, B S (X, X*) contains no operators that 
multiplied by a large enough number become larger than any given operator, like the identity 
/ £ B S (H). One can show that any such operator E must satisfy (Ex,x) > a||x || 2 for some a > 0, 
and therefore (Ex,x) 2 is an equivalent Hilbert norm on A". 

Moreover, Bf(X,X*) may not be generating, i.e. not every bounded and symmetric operator 
can be represented as a difference of two positive definite operators. The authors of [45] prove that 
Bf(X, X*) is generating if and only if every P £ B S (X,X*) factors through a Hilbert space (by 
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the reproducing kernel Hilbert space construction any P G Bf(X,X*) always so factors, see [131 
2.1] and Section [2]). An example of non-factorizable P G B s (l p , I*) for 1 < p < 2 is given in 
hence Bf(l p ,l*) is not generating for such p. Moreover, any infinite-dimensional L p (p) contains 
a complemented copy of l P EH Thm.3.3], so Bf(L p (p), L p (y)*) is not generating for 1 < p < 2 
either. However, some properties of the order are shared with B S (H), as the next theorem shows. 

Theorem 4. Let || • || denote the induced operator norm on B(X, X*), then 

(i) If P G B S (X,X*) then ||P|| = sup \(Px,x)\, hence || • || is monotone on Bf(X,X*); 

\\x\\=l 

(ii) IfPe B+(X,X*) then \\Px\\ 2 < ||P|| (Px,x) for all x G X; 

(iii) order bounded monotone sequences in B S (X,X*) strong operator converge. 

Proof, (i) The proof is analogous to the Hilbert space case, see e.g. [27, Ch.VI], Let M : = 
sup || 3 ,|| =1 \ (Px,x)\, obviously M < ||P||. By symmetricity of P and the parallelogram identity, 

4 (Px,y) = ( P(x + y),x + y) - (P(x - y),x - y) < \(P(x + y),x + y)\ + \(P(x - y), x - y)\ . 

Therefore, 

4 (Px,y) < M (\\x + y\\ 2 + ||x - y\\ 2 ) < 2 M (||a;|| 2 + \\y\\ 2 ) . 

Cancel 2 and set y = jj^Px in the last inequality to obtain 

2 ( Px 'wk Px ) = nPA 11x11 - M ( IW|2 + H Pl H 2 )= 2 M INI 2 ■ 

Canceling 2 ||a;|| we get ||Pr|| < M ||x|| for all x, and hence the converse inequality ||P|| < M 
holds. 

Recall that a norm is called monotone if 0 < P < R implies ||P|| < ||P||. But if P < R then 
(Px, x) < ( Rx , x) for all x G X by definition of order. For positive definite P the absolute value 
in the formula just proved can be omitted, so taking supremum over all x with ||x|| = 1 gives 

\\P\\ < |» 

(ii) This inequality is derived in passing in m- Consider tp G X*, for any £ > 0 there is a 
normalized element g G X such that (<p,g) > (1 — e:)||<^|| since ||</?|| = sup {p,x). Setting 

:eEX,||:e||=1 

g := \\<p\\g we have ||<?|| = ||(^|| and (<p,g) > (1 — e)||(/?|| 2 , in other words g approximates in X the 
tangent functional to <p (which in general belongs to X**). 

Applying the last inequality to <p = Px we have (1 — £)||Pa :|| 2 < ( Px,g). Moreover, if P > 0 
then (Px, y) is a semi-definite inner product on X, so it satisfies the Cauchy-Schwarz inequality 
\(Px,y)\ 2 < (Px,x) ( Py,y). Setting y — g and combining the inequalities we obtain 

(l- e) 2 \\Px\\ A <\(Px,g)\ 2 < (Px,x)(Pg,g) 

<||P|| |M| 2 <P*,z) = ||P|| \\Px\\ 2 (Px,x). 
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Canceling ||Px|| 2 we get (1 — £) 2 ||Px|| 2 < ||P|| (Px,x) (if Px = 0 this inequality holds trivially), 
and since e > 0 is arbitrary the desired inequality follows. 

(iii) For any x G X the sequence (P n x, x) is monotone and bounded, and hence converges by the 
Weierstrass theorem. By the polarization identity a limit l x (y) of (P n x,y) exists for any x, y e X. 
Multiplying P n by —1 and adding a lower bound to all of them if necessary we may assume that 
P n are positive definite and monotone increasing. Since by (i) the operator norm is monotone the 
norms of P n are uniformly bounded, ||P n || < M. This implies that l x is a continuous functional on 
A", and one can define a linear operator from X to X* by setting Px := l x . 

By construction, P n weak operator* converge to P. Moreover, by (ii) 

||(P„-P)x|| 2 < ||P„-P|| ((P n — P)x, x) < 2M((P n — P)x, x) -> 0 . 

n—»• oo 

Thus, P n strong operator converge to P. □ 

Operators in B S (X,X*) can not be composed directly, but any operator N 6 B s {X*,X) can 
mediate composition and define a product on B S (X, X*) by (P, R) H y PNR. Since \\(P n NR n — 
PNR)x || < ||P n || || N\\ || (R n — P)x|| + || (P„ — P)NRx\\ this product is sequentially strong opera¬ 
tor continuous. Therefore, we have the following corollary for the nonlinear part of the Riccati 
equations. 

Corollary 2. For any N € B s (X*,X) the product on B S (X,X*) defined by (P, R) i-G PNR is 
continuous on order bounded monotone sequences. 

4 Concave equations and stabilizing solutions 

We now have the analytic tools to construct solutions to the Riccati equations. However, it turns 
out that the construction mostly uses order properties, and can be abstracted from the Riccati 
specifics. What we need is that the equations are quasi-linear with the linear part generating a 
positive semigroup, and with the non-linear part being concave. This is the generality adopted in 
this section. 

We start by recalling some basic facts about ordered Banach spaces, positive semigroups 
on them, and the order concavity. We also recall some key results of [19] on abstract Lyapunov 
equations. They are then used to construct a special class of solutions to quasi-linear concave 
equations by Newton’s iteration. In dimension one Newton’s iteration provides monotone conver¬ 
gence of approximations for concave equations not just locally, but globally. We prove a similar 
result in ordered Banach spaces. 

Let A be a Banach space with predual X * partially ordered by a closed proper cone X + , see [6], 
App.2], The predual cone is defined in the usual way Xfi := e X* \ (ip, x) > 0 for all x G X + }. 
Recall that a Cq semigroup T{t) on A is a weak* continuous semigroup adjoint to a C 0 semigroup 
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7*(t) on A’*. Its generator A is the weak* limit of the difference quotients, which is weak* densely 
defined, weak* closed, and adjoint to the generator A* of the predual semigroup T(t) |5], 3.1.2], If 
Tit) is positive, i.e. preserves X + , then so is T(t), and vice versa. 

We will be studying quasi-linear equations of the form Ax + $(x) = 0 on X + , where A is 
the generator of a positive Cq semigroup, and $ is Gateaux differentiable [16:, ch.XVII. 1.1],[32l 
1.2.1], and concave in some sense. Algebraic operator Riccati equations are not the only ones of 
this form, so are some reaction-diffusion equations in mathematical chemistry for example. With 
the help of Newton’s iteration the problem reduces to solving a sequence of linear equations of 
the form Ax + z = 0 . Our assumptions about A and z turn them into abstract analogs of the 
operator Lyapunov equations. The main technical tool is a criterion that relates solvability of 
Ax = —z on X + to positive invertibility of —A, and to exponential stability of T(t), if z satisfies 
a non-degeneracy condition below. 

Definition 7. An element z G X + is called an L 1 detector for T(t) on a subset T C Xff if for 
every p G T: 

/*oo POO 

/ {T(t)p, z) dt < oo => / \\T(t)p\\ dt < oo . 

Jo Jo 

A subset T C Xff is called a stability subset for a class of semigroups if for every semigroup in 
the class J 0 °° ||71(£)v?|| dt < oo for all p G T implies that the same inequality holds for all p G Xff. 

Basically, the detector condition asks that z correctly ’’detect” asymptotic behavior of Tit) on 
every element of T. The terminology comes from examples in control theory, see Section [5l Readers 
familiar with order unit spaces [ 6 ], App. 2 ] will immediately see that any order unit is an L 1 detector 
on all of Xjf for any semigroup. But semigroups can have L 1 detectors even in spaces without 
order units. It is also immediate from the definition that if w > z and z is an L 1 detector, then 
so is w. Note that L 1 stability of T{t) on T does not in general imply L 1 stability on the entire 
*. + , even if T is dense in it. Of course, T = Xff is a stability subset for all Cq semigroups, but 
the class of Lyapunov semigroups admits a much smaller stability subset, see Lemma [3l 

We will assume that Xff is generating, i.e. every element in X * is a difference of two positive 
ones. We do not however impose this condition on X + keeping in mind our example of X = 
B S {X, X*) and X * = X®^X. If J 0 °° ||7^(i)<£>|| dt < oo holds for every p G Xf~ and Xf~ is generating, 
then it holds for every p G X *. Then by the Datko-Pazy theorem |4J II.1.2.2], [ 6 j Prop.9.4] we have 
||7^(t)|| = ||T(f)|| < Me~ £t for M ,e > 0 since Tit) is a C 0 semigroup, i.e. T(t) is exponentially 
stable. This makes the following result less surprising. We state it for Cq semigroups, which is all 
we need here, but in [191, Thm.3] it is proved in greater generality. 

Theorem 5. Let T(t) be a positive Cq semigroup on X with the generator A and an L 1 detector 
z on a stability subset T. If Xff is generating the following conditions are equivalent: 

(i) Ax = —z has a positive solution x G D_ 4 D X + ; 
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(ii) Pit) is exponentially stable; 

(iii) A has a bounded inverse on X and —A^ 1 > 0. 

This is an abstract version of Wonham’s stability criterion for Lyapunov equations. The most 
non-trivial implication is from (i) to (ii). When using the theorem the hardest part to verify is the 
detector condition. One way is to use continuous final observability that implies L 1 detectability 
m Thm.2], and can often be established in applications by priori estimates, see e.g. [7J 3.27]. 

For matrices there is a notion of matrix concavity/convexity [3], and it was used implicitly 
by Klcinman in the original proof of monotone convergence of approximations na to solutions to 
algebraic matrix Riccati equations. Its use is made explicit for a larger class of matrix equations in 
| 8 ], and Zabczyk generalized Kleinman’s implicit approach to operator Riccati equations in Hilbert 
spaces [35]. Ordinary and matrix concavity are particular cases of order concavity defined next. 

Definition 8. Let X be an ordered Banach space. A map <f> : X I) D —> X is called order concave 
if for all a G [0,1] and x,y E X + , we have <F(ax + (1 — a)y ) > a$(x) + (1 — a)$(y). 

If D is a convex set (in the usual sense) and $ is Gateaux differentiable one can prove as in the one 
dimensional case that order concavity is equivalent to <F(y) — <F(x) < &'(x)(y — x) for all x,y G D, 
a condition which is often simpler to verify. The Riccati map < h(T > ) = —PNP + Q is Gateaux 
differentiable, and when IV > 0, order concave on B S (X, X*), see Section [5j 

In the rest of this section A denotes the generator of a positive Cq semigroup on X with 
the domain Z/ 4 , and $ : X + —> X is an order concave Gateaux differentiable map. We set 
F(x) : = Ax + $(x) for x G an d for all x, y G X + we denote for convenience: 

F\x) :=^ + $ / (x) 

L(x) := <L(x) — $'(x)x (8) 

T(x, y) := <L(x) — <f>(y) + &(x)(y — x). 

One can see that for x G T .4 changing $ to F in the definitions of L and T does not change 
the result. Here T measures deviation of $ from its tangent, and in dimension one L(x) = 
$(0) + T(x,0) is closely related to the Legendre transform of F and $. The order concavity is 
equivalent to T(x, y) > 0 for all x, y G D, so if 0 G D then L(x) > $(0) for all x G D. 

Newton’s iteration is a classical method for solving nonlinear equations of the form F(x) = 0, 
see e.g. PI ch.XVIII]. By construction, the iteration equation is F(x n ) + F’(x n )(x n+ 1 — x n ) = 0. 
It can be rewritten as F\x n )x n+ 1 = —L(x n+ 1 ), which is of the form Ax = —z. The next lemma 
describes the effect of performing one step of Newton’s iteration in our setting. 

Lemma 2 (Iteration step). Assume that X* is generating and F'(u ) \= A + Q'(u) generates a 
positive Cq semigroup TF'{ u )if) for every u G X + . Suppose $(0) > 0 and for some x G X + the 
semigroup TF'{ x )(t ) is exponentially stable. Then: 
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(i) There is a unique solution y E T>a D X + to F'(x)y = —L(x); 

(ii) F(y) = -V(x,y) < 0,- 

(iii) U moreover y E T>a D X + and F(x) < 0, then y < x; 

(iv) If moreover L(y) is an L 1 detector for TF'( y ){t) then TF'( y )(t ) is exponentially stable. 

Proof, (i) By assumption we have L(x) > $(0) > 0. Since F'(x) is stable —F'(a ;) _1 > 0 by 
Theorem 0 and y = — F'(x)~ 1 L(x) E T>f'(x) C X + = T>a H X + . 

(ii) By definition of F and the equation for y we compute F(y) = F'(x)y + L(x) — \l t(x,y) = 
— T(x, y) < 0. 

(iii) Since x E V>a the equation for y is equivalent to F(x) + F'(x){y — x) — 0. And since 
—F'(x ) _1 > 0 by TheoremO y — x = —F'(x)~ 1 F(x) < 0. 

(iv) Setting y — x in the identity from (ii) we obtain F'(y)y = —(L(y) — F(y)). This is an 

abstract Lyapunov equation to which y is a solution. Since —F(y) > 0 by (ii) and L(y) is an 
L 1 detector L(y) — F(y) is also an L 1 detector. Since y > 0 our equation has a positive solution 
fulfilling part (i) of Theorem [5] Part (iii) of that theorem gives the desired conclusion. □ 

One can see that in conditions of Lemma [2] the assumptions imposed on x are reproduced and 
even improved for the next approximation y. Therefore, assuming a suitable initial guess Xq can 
be found, we can generate a monotone decreasing positive sequence x n that (hopefully) converges 
to a solution to F(x) = 0. One kind of convergence already follows from our assumptions. If A* 
is generating then monotone bounded sequences always weak* converge in X. Indeed, numerical 
sequences (ip,x n ) converge by the Weierstrass theorem for any ip E Xf~, and therefore for any 
ip E A*. So x n is weak* Cauchy. Since bounded sets in the dual to a Banach space are weak* 
precompact by the Alaoglu theorem we conclude that x n weak* converge to some Xoo E X + . 

However, weak* topology is so weak that nonlinear maps of interest are rarely continuous in 
it. Fortunately, monotone sequences can often be proved to converge in a stronger sense, in which 
the relevant maps are continuous, see Corollary [2l 

Definition 9. A map 0 : X D D —> X is called continuous on bounded monotone sequences if 
for any bounded monotone sequence x n with weak* limit x the values 0 (x n ) converge to 0 (x) as 
n —» oo. 

To simplify the statement of the following theorem it is convenient to introduce some termi¬ 
nology. The motivation comes from control theory and is explained in Section 0 

Definition 10. We say that an element x E X + stabilizes {A, $) if the operator F' (x) := A+$'{x) 
generates an exponentially stable positive Cq semigroup. A pair (^4, <3>) is called exponentially 
stabilizable if some x E X + stabilizes it. A pair ( A , <L) is called globally L 1 detectable if for every 
x 6 X + the element L(x) := $(x) — d >, (.T)a: is an L 1 detector for Tf'( x ){I) on a stability subset T 
(we suppress the dependence on T in the name). 
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The next theorem is an abstract version of our main result. 

Theorem 6. Let a pair (A , <3>) be exponentially stabilizable and globally L 1 detectable, <h(0) > 0 
and x $>(x), (x,y) i->- $'{x)y be continuous on order bounded monotone sequences. Pick xq G 
X + so that T F g X0 )(t) is exponentially stable. Then the sequence defined recursively by F(x n ) + 
F\x n )(x n+ 1 — x n ) = 0 is well-defined, belongs to T>aFX + , and is monotone decreasing for n > 1 . 
It weak* converges to a unique positive solution to F(x) = 0, and this solution is stabilizing. 

Proof. Existence: By solving the iteration equation and applying Lemma [2] we produce x n G 
T>a H X + with xi > x 2 > • • • > 0. The iteration equation can be rewritten as Ax n+ i = —$(a; n ) + 
&(x n )x n — &'(x n )x n+ i. Since x n is a bounded monotone sequence there is a weak* limit x n -» 

n—¥oo 

Xoo G X + . By continuity assumptions on <L: 

Ax n+l -* -$(iCoo) + ^'Ooo^oo - ^'Ooo^oo = -$(^oo)- 

n —>00 

Since A generates a Cq semigroup it is weak* closed, so Xoo G Va and Ax 00 = —^(xoo). Thus, 
F(x 0 o) = Ax oo + ^(xoo) = 0, and x ^ is a solution. 

Stabilization: Let x G T>a fl X + be any positive solution to F(x) = 0, then F\x)x = 
— (0 — F'(x)x) = —(F(x) — F'(x)x) = —L(x). Since L{x) is an L 1 detector for F\x) by global 
detectability, we conclude from Theorem [5] that TF<( Xao )(t) is exponentially stable. Thus x, and x^ 
in particular, is stabilizing. 

Uniqueness: For any positive solution x we have by concavity 0 = F(x) — F(xoo) < F , (x 00 )(x — 
Xoo). Since TFq Xoo )(t ) is positive and exponentially stable —F'(x oo) -1 > 0. Therefore, x — Xoo < 
— F , (x oo ) _1 0 = 0. Switching the roles of x and x^ we also get Xoo — x < 0, and x = Xoo ■ □ 

The condition F( 0) = $(0) > 0 can be somewhat relaxed to F(x) > 0 having a solution 9 G X + 
and a stabilizing x 0 existing with x 0 > 9. This case reduces to ours by passing to F{x) : = 
A{x + 9) + <L(a: + 9), which satisfies the conditions of Theorem [6l If Xqo is a solution to F(x) = 0 
then Xoo ■= Xoo + 9 is a solution to the original F(x) = 0. In this case Xqq is still a unique stabilizing 
positive solution, but other positive solutions may exist. They are smaller than or incomparable 
to 9. Since any of them can be chosen in place of 9 we conclude that Xqo is also the largest positive 
solution to F(x) = 0, and even to F(x) > 0. Such characterization of the stabilizing solution is 
also valid for stochastic matrix Riccati equations p]. 

5 Stabilizability and detectability 

In this section we prove our main result by expressing the conditions of the abstract Theorem [6] 
explicitly for the Lyapunov semigroups, and then verifying them. We also recall stabilizability and 
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detectability conditions from control theory, and explain their relation to the abstract versions 
from the previous section. 

One of equivalent ways to define detectability [33; 3.6] in finite dimensional spaces is to call 

—y 0. In infinite 


UA 


x 


a matrix pair (C, A) detectable if for every x: Ce tA x - > 0 implies e' 

t—t oo t —^OO 

dimensions different types of convergence to 0 are no longer equivalent, so many generalizations 
are possible, ours is one of them, cf. [3] II.1.2.2], 

Definition 11. Let X be a Banach space and T(t ) be a Co semigroup on it with the generator 
A. Let V be a Hilbert space and C : X —» V be a bounded operator. The pair ( C , A) is called 
detectable in L 2 if for all x G X 


||CT(i)a ;|| 2 dt < oo 


||T(f)a ;|| 2 dt < oo . 


Recall that any Lyapunov semigroup is adjoint to 7 *(t)(x <g> x) = T(t)x <g) T[t)x. It follows 
that ||T(i)|| = 11 %(t) || = ||T (£)|| 2 and exponential stability of T(t) is equivalent to that of T(t). A 
crucial observation is that (X<g) 7r A) + has a very simple stability subset for Lyapunov semigroups 
as a consequence of the Datko-Pazy theorem [3], II. 1.2.2], [ 6 ], Prop.9.4], The theorem states that 
for Cq semigroups if J 0 °° \\T(t)x\\ p dt < oo for some p > 1 and all x G X then T(t) is exponentially 
stable. 

Lemma 3. The set of tensor squares T := {x <g) x| x G X} is a %(t) invariant stability subset 
of (A"(g) 7r X) + for the class of Lyapunov semigroups. Moreover, Jo°° l|7*(^)p|| dt < oo for all p G T 
implies even exponential stability of 7l(t). 

Proof. The invariance is obvious from the definition of 7*(t). And if / 0 °° ||7^(t)(x <g> x)\\dt = 
/ 0 °° \\T(t)x®T(t)x\\ dt = / 0 °° ||T(t)x || 2 dt < oo for all x G X then by the Datko-Pazy theorem T(t), 
and hence 7 *(t), is exponentially stable. □ 

Setting R := C*C we can rewrite the definition of L 2 detectability as 


((7 ~*(t)(x <g) x), R )) dt < oo 


|| %{t){x <g) x)\\ 2 dt < oo . 


(9) 


In other words, L 2 detectability of (C , A) is equivalent to R being an L 1 detector in the abstract 
sense of Definition [7] for the Lyapunov semigroup of T{t) on the stability subset of tensor squares. 
We will also say that R is an L 2 detector for the generator A of T(t). 

But for (Ca, $) to be globally L 1 detectable in the sense of Definition flOl we will need L(P) = 
Q + PNP to be an L 1 detector of Ca-np f° r every P G Bf(X,X*). Simple L 2 detectability 
of ( C, A) is not enough to ensure detectability globally. One can show however that exponential 
detectability of ( C , A) is enough. A better alternative would be to prove that L 2 detectability is 
preserved by the Newton’s iteration, but this seems unlikely to hold in general. 
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Lemma 4. Let A be the generator of a positive Co semigroup T(t) on X, B G B(U,X) and 
C G B(X,V) for some Hilbert spaces U, V. Suppose that (C, A) is exponentially detectable. Then 
for any I\ € £>(W, [/) £/je operator C*C + K**K is an L 2 detector for A — BI\. 

The proof is analogous to the Hilbert case [35) Lem.3], see also [U Thm.I.1.2.6], and is omitted. 

Corollary 3. Suppose A is the generator of a Co semigroup on X, and $(P) —PBB*P + C*C 
for B G B(U, X), C G B(X, V) with some Hilbert spaces U, V. If ( C , A) is exponentially detectable 
then (£^ 4 ,$) is globally L 1 detectable on Bf (X, X*). 

Proof. By Definition ITOl we need to show that for every P G Bf(X,X*) the element L(P) = 
$(P) — &(P)P = PBB*P + C*C is an L 1 detector for the Lyapunov semigroup generated by 
Ca + $ 7 (P) = jCa-bb*p, see (fill) and Corollary [H Setting K := B*P G B(X, U ) and noting that 
K*K = PBB*P is symmetric we see that this is exactly the conclusion of Lemma 01 □ 

Setting K = 0 in Lemma H] one can also see that exponential detectability implies detectability 
in L 2 . Specializing Theorem [5] to our current setting we prove Theorem [Tj from Section [2] next. 
It is a generalization of the stability criterion in terms of Lyapunov equations proved for reflexive 
spaces in [T9], the matrix case is due to Wonliam [33) 12.4]. In fact, we prove a stronger claim that 
only asurnes that (C, A) is detectable in L 2 . 

Proof of Theorem [ 1 } By inspection and Lemma [3] X = B S (X, X*), T(t), z = C*C and T = 
(X<g)X) + satisfy the conditions of Theorem [5] We see that Theorem [5](i) is equivalent to the (i) 
here. Claims (ii),(iii) of Theorems [5] are equivalent to (ii),(iii) here as well since the exponential 
stability of T(t ) is equivalent to that of T(t). □ 

We now turn to stabilizability. The map <I>(P) = —PNP + Q is clearly Gateaux differentiable 
with T'fP)/? = — PNP — RNP , and the auxiliary maps from (| 8 ]) become 

L(P) = 4>(P) - 4>'(P)P = PNP + Q (10) 

T(P, Q) := $(P) - $(R) + &(P)(R - P) = (P - R)N(P - R ). 

The last formula implies in particular that <3? is order concave when N > 0 since T (P, Q) > 0 is 
equivalent to the definition of order concavity for Gateaux differentiable maps. For the full Riccati 
map F(P) = CaP + $(P) the formal derivative is equal to 

F'(P)R = C a R + &(P)R = A*R + RA — PNR - RNP = (A* - PN)R + R(A - NP). 

Although P* yl P in non-reflexive spaces we still have ( B*P)* = PB because the domain of B 
is a Hilbert space, so B** = P, and Ran (P) C X while P*\ x = P. Therefore also (NP)* = 
(BB*P)* = (B*P)*B* = PBB* = PN, and using Corollary [2] we conclude that 

F'(P)R = (A - NP)*R + R(A - NP) = jC A -npR ■ ( 11 ) 
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This means that (Ca, $) is exponentially stabilizable in the sense of Definition [TO] if and only if 
there is P G Bf(X,X*) such that A — NP is exponentially stable. For N = BB* we will reduce 
this to existence of K G B(X, X *) such that A — BK is exponentially stable, which is the standard 
definition of a pair (A, B) being exponentially stabilizable [33] 2.4], [36] 1.2.5]. 

Lemma 5. Let A be the generator of a positive Co semigroup on X, B G B(U,X ) and C G 
B(X, V ) for some Hilbert spaces U. V. Suppose that (C , A) is exponentially detectable and set 
<f>(P) := —PBB*P + C*C. Then ( Ca , 4>) is exponentially stabilizable in Bf (X. X*) if and only if 
(A, B) is exponentially stabilizable. 

Proof. In this proof the lower indices under both semigroup symbols T(t) and T (t) indicate the cor¬ 
responding generator of T(t). Since Ca + &'(P) = Ca~bb*p one direction is trivial. If Ta- BB*p(t ) 
is exponentially stable then so is Ta~bb*p(1 ) and one can set K := B*P. 

Conversely, let Ta-bk^), and hence Ta-bk( t), be stable for some K G B(X,U). Consider 
the Lyapunov equation Ca-bkP = —(C*C + K*K ). The cone (W<g) vr W) + is generating and by 
Lemma 0] C*C + K*K is an L 1 detector for Ta-bk( t), so by Theorem [5] there exists a solution 
P G Bf(X,X*). By Corollary [2] we have 

Ca-bb*pP = Ca-bkP + Cbk-bb*pP 

= ~(C*C + K*K) + (BK - BB*P)*P + P(BI< - BB*P ) 

= -C*C - K*K + K*B*P - PBB*P + PBI< - PBB*P 

= ~(C*C + PBB*P+ (K - B*P)\K - B*P) s j. (12) 

This is another Lyapunov equation that P solves. By Corollary [3l C*C + PBB*P is an L 1 detector 
for Ca~bb*p, all the more so after adding (K — B*P)*(K — B*P ) > 0. Existence of a solution 
P > 0 guarantees by Theorem [5] that Ta-bb*p(P) is exponentially stable, so this P provides the 
required stabilization of (Ca, $)• □ 

We are ready to prove the main theorem. 

Proof of Theorem [2] Let Ca be the Lyapunov generator corresponding to A, i.e. C \P = A*P + 
PA on T>a, and let < F(P) := —PBB*P + C*C. By Corollary [3] and Lemma O the pair (Ca, 4>) is 
exponentially stabilizable and globally L 1 detectable. The choice of P 0 is justified by Lemma 0 
Obviously, <F(0) > 0, and P ^ $(P) and (P, R ) t-)- <F'(P)P = — PBB*R — RBB*P are continuous 
on bounded monotone sequences by Corollary[2] Therefore, by Theorem [6] there is a unique positive 
definite solution P to the Riccati equation, which is stabilizing and a weak* monotone limit of 
solutions P n to (J6]) . Since P n is monotone decreasing and bounded by Theorem [4](iii) it also strong 
operator converges to the same limit. It remains to prove the quadratic convergence estimate. 
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Let P(P) := A*P + PA — PBB*P + C*C on T>a■ By Lemma [2](ii) we have P(P n+ 1 ) = 
—T(P n ,P n+ 1 ); where T is as in ([8]) with N = BB*. Therefore, by concavity and (JHI) 

F\P){P n+1 - P) > P(P n+1 ) - P(P) = -T(P n , P n+1 ). 

Since P'(P) = Ca~bb*p and P is stabilizing the semigroup generated by Ca~bb*p is exponentially 
stable, and the inverse — F'(P)^ 1 is bounded and positive by Theorem [fDpii). Therefore, 0 < 
P n _ |_i — P < — F , (P) _1 T(P n , P n+ 1 ). Since the norm on B S (X , X*) is monotone by Theorem |4^i) we 
estimate 

l|F„ + i-P|| < ||F'(F)- 1 ||||>l<(F„,P„ + i)||. 

It follows from (18]h the inequalities 0 < P n — P n +\ < P n — P, and the monotonicity of the norm 
that 


||T(P n ,P n+1 )|| = ||(P n -P n+ 1 )T?£r(P„~P ri+1 )|| < ||P|| 2 ||P n -P n+1 || 2 < ||P|| 2 ||P n -P|| 2 . 

Setting k := ||(£ j 4 _ bb *p) -1 ||||P|| 2 we obtain the estimate. □ 
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